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t — ' Using the Arcminute Microkclvin Imager (AMI) at 16 GHz and the Very Small 

\ Array (VSA) at 33 GHz to make follow-up observations of sources in the New Ex- 

0^ . tragalactic WMAP Point Source Catalogue, we have investigated the flux-density 

variability in a complete sample of 97 sources over timescales of a few months to 
* 1.5 years. 

We find that 53 per cent of the 93 sources, for which we have multiple observations, 
y\ ' are variable, at the 99 per cent confidence level, above the flux density calibration 

^ , uncertainties of « 4 per cent at 16 GHz; the fraction of sources having varied by 

more than 20 per cent is 15 per cent at 16 GHz and 20 per cent at 33 GHz. Not 
only is this common occurrence of variability at high frequency of interest for source 
physics, but strategies for coping with source contamination in CMB work must take 
this variability into account. 

There is no strong evidence of a correlation between variability and flux density 
for the sample as a whole. For those sources classified as variable, the mean fractional 
r.m.s. variation in flux density increases significantly with the length of time separating 
observation pairs. Using a maximum-likelihood method, we calculate the correlation 
in the variability at the two frequencies in a subset of sources classified as variable 
from both the AMI and VSA data and find the Pearson product-moment correlation 
coefficient to be very high (0.955 ± 0.034). We also find the degree of variability at 
16 GHz (0.202 ± 0.028) to be very similar to that at 33 GHz (0.224 ± 0.039). 

Finally, we have investigated the relationship between variability and spectral 
index, af| g 5 (where S oc ^~ Q ), and find a significant difference in the spectral indices 
of the variable sources (—0.06 ± 0.05) and non- variable sources (0.13 ± 0.04). 
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1 INTRODUCTION 

The existence of the time variability of the flux densities of 
some extragalactic radio sour ces has been established since 
the 1960s (see e.g. Dent 1965). Long-standing issues, stem- 
ming mostly from observations at < 5 GHz, have concerned 
the effects of the path between the source and us - with im- 
plications for the propagation medium - and of the possible 
roles of variability of the source luminosity - with implica- 
tions for how the sources work. More recently, as studies of 
structure in the CMB have become possible (requiring ob- 
servations at > 10 GHz), knowledge of variability at high 
frequency is becoming essential in order to remove CMB 
foreground contamination. We provide some illustration of 
the importance and causes of variability, as follows. 



1.0.1 Interstellar scintillation 



IShapirovskaval (1 19781 ) and lRickett. Coles fc Bourgoisl (|l984l ) 
suggested that refractive interstellar scintillation (ISS) could 
be responsible for the intra-day variability of extragalactic 
radio sources at v < 5 GHz. The basis of the mechanism 
is that the refractive index, n, from an element of plasma 
in a cloud in the interstellar medium (ISM) in our galaxy 
has the usual form n = [1 — ( — ) 2 ] 1 ^ 2 , where u is the obser- 
vation frequency and v v is the plasma frequency, but that 
the variation in n in the plane of the sky across the cloud 
(or indeed resulting from more than one cloud in the line of 
sight) produces scintillation as the galactic ISM moves with 
respect to both us and the radio source. The form of n(y) 
indicates that these scintillations will not be important at 

high enough v. 

Many studies (see e.g . Gregorini. Ficarra fc Padriellil 



ll986l ; ICawthorne fc Ricket^l 19851 ; Tspanglei • et allll989l ) sup- 
port this notion in that they show that variability tends to 
increase as the galactic latitude \b\ of t h e sour ces decreases. 
Further, iDennett-Thorpe fc de Bruvnl (|2002|l following a 
techni que due to lLang fc Rickettl 1 1970T ) and Jauncev et all 
(2000), used two telescopes at different locations on Earth 
to demonstrate the effect on the flux density at 4.9 GHz 
of J1819+3845 of the Earth moving through a scintilla- 
tion pattern at the Earth caused by the local ISM. Finally, 
iLinskv. Rickett fc Redfieidl (2008) have shown unequivocally 
that the motion of individual clouds in the local ISM, espe- 
cially as cloud edges pass the line of sight, is an important 
mechanism of variability at 5 GHz. 



1.0.2 Intrinsic source variability 

It is also clear that there is a component of variability that 
is intrinsic to some sources. The existence of this type of 
variability is fully supported by VLBI observations of the 
motion and lumino sity evolution of synchr otron emitting 
components (see e.g. lZensus fc Pe"a rson 1987). For a moving 
source, S t>, the observed flux density, is related to S, the 
flux density that would be observed at the same frequency 
v in the comoving frame, by S \>[y) = S(u)D i+a , where a 
is the spectral index and D the boost factor which depends 
on both the speed of the source and its angle to the line of 



sight. (This equation is true for optically thin sources and 
spherical optically thick sources.) For a relativistically mov- 
ing source, curvature in the tra jectory can therefore prod uce 
flux density variation (see e .g. Blandford fc Konigllll979T) . 

iRickett. Lazio fc Ghigol l|2006T ) have discussed the de- 
gree of variability as a function of v. They studied 146 ex- 
tragalactic compact radio sources (a very important but not 
a complete sample of all sources above a given flux density) 
monitored over 1979-1996 with the Green Bank Interfer- 
ometer at 2 and 8 GHz; major conclusions are that the ISS 
component of variability is much stronger at 2 GHz than 
8 GHz, and that intrinsic variability is of prime importance 
at high radio frequency. 

We now focus on variability at > 15 GHz. Given the 
foregoing work, we expect this variability to be dominated 
by intrinsic luminosity variability - so at 15 GHz we are 
dealing with the intr insic properties of the radio sources - 
although the work of ILinskv. Rickett fc Redfieidl clearly im- 
plies that a small fraction of sources will have substantial 
variability caused by the ISS as the edges of ISM clouds 
pass the line of sight between us and the radio source, even 
at high v. 

1.0.3 Our work 

There are two large-scale surveys of extragalactic radio 
sources that are complete in flux density at > 15 GHz and 
which study variability: the 9th Cambridge (9C) su rvey com- 
plete to » 25mJy at 15 GHz jBolton et al.ll2006t ) and the 
Australia Telescope Compact Ar ray (ATCA) 20-GH z pilot 
survey complete to « 100 mjy (Sadl er et all 120061) . Here, 
we study variability in a high- flux-density (> 1 Jy), very- 
large-area survey. As well as being complementary to the 
two other studies, our work has implications for forthcom- 
ing CMB work. 

Our work follo ws on from 

I AMI Consortium: Davies et all (120091 ') (Paper I) which 
employs the New Extragalactic WMAP Point Source 
(NEWPS) catalogue and is complete to 1.1 Jy at 
33 GHz. Using the Arcminute Micr okelvin Imager (AMI; 
lAMI Consortium: Zwart et all |200S|) at 16 G Hz and the 
Very Small Array (VSA: IWatson etail l2003h at 33 GHz, 
we have carried out simultaneous observations of a subset 
of the sample to investigate variability over timescales 
ranging from a few months to w 1.5 years. All the AMI 
data presented in this paper were obtained using the AMI 
Small Array. Observation and reduction strategies are as 
described in Paper I. 



2 FLUX DENSITY VARIABILITY 

2.1 Variability of the source population as a whole 

2.1.1 Previous work 

iBolton et all |2006) studied the variability of 51 faint 
sources from the 9C survey (complete to w 25 mjy at 
15 GHz) at 15 GHz on a 3-year timescale. They defined the 
fractional variation as 
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Table 1. Distribution of the unbiased varia bility index V TI 
20 GHz for sources in the lSadler et al . (2006) sample. 



at 



Vims (per cent) 


n 


Fraction (per cent) 


< 10 


63 


58 ±7 


10-20 


29 


27 ±5 


20-30 


11 


10 ±3 


> 30 


5 


5±2 


Total 


108 





where Xi are individual flux density measurements for the 
same source, n is the number of data points and \i is the 
mean flux density. For the 17 sources in lBolton et al.l classi- 
fied as variable, the fractional variation ranges between 8.4 
and 70 per cent, w ith a median of 14 per cent. 

ISadler et all l|2006j ) studied the variability of 108 
sources from the ATCA 20-GHz pilot survey (complete to 
100 mjy at 20 GHz) at 20 GHz on a 1-2 year timescale. For 
this purpose, they used follow-up observations at two epochs 
during 2003-04. The y quantified the variability in the same 
way as lBolton et all except that they took the uncertainties 
in individual flux density measurements into account. They 
defined the unbiased variability index, Kms by 



(2) 



where o~i are errors on indiv idual flux density measure - 
ments. This expression follows lAkritas fc Bershadvl (l996), 
who make it clear that the value of <7j must be the same 
for each of the obser vations of a particular source. For all 
sources, Sadler et al.l obtained a median value of Kms at 
20 GHz of 6.9 per cent. Table [1] shows their distribution of 

Vrms • 



2.1.2 Methods 

Good knowledge of the uncertainties in flux density mea- 
surements is crucial for the study of variability. Great care 
was put into determining these before taking them into ac- 
count in quantifying the variability (see Paper I). Because 
the VSA is not equatorially mounted, the linear polarisation 
measured with the VSA changes as a function of parallactic 
angle. However, we almost entirely overcome this difficulty 
by observing the sources at very similar HAs in different ob- 
serving runs. We also note that, if polarisation were playing 
an important role, we would have expected to see greater 
scatter in the measured flux densities for each source within 
a single observing run. Moreover, in general the WMAP- 
detected point sources are not strongly pola rised: for the 
source s detected in the WMAP 5-year maps, IWright et al.l 
(2009) have measured a mean polarisation at 33 GHz of just 
2.2 per cent, although we have not determined whether there 
is a relation between polarisation and variability. We there- 
fore do not expect polarisation to be an important factor in 
the variability observed with the VSA. 

We have tried to analyse the variability in such a way 
as to make our re sults comparable to those obtained by 
Sadler et al.l l|2006h . So, we have only made use of the first 
and last flux density measurements from both telescopes for 
each source, yielding similar time intervals. In the AMI data, 



Table 2. Distribution of the unbiased variability index V at 16 
and 33 GHz for sources in our sample. The errors on the fractions 
are Poisson errors. 



V (per cent) 




n 


Fraction (per cent) 




16 GHz 


33 GHz 


16 GHz 


33 GHz 


< 10 


61 


55 


66 ± 8 


59 ± 8 


10-20 


18 


19 


19 ±5 


20 ±5 


20-30 


8 


12 


9 ± 3 


13 ±4 


> 30 


6 


7 


7±3 


8±3 


Total 


93 


93 
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Figure 1. Unbiased variability index V at 16 against that at 
33 GHz. 



time intervals range between 1 and 18 months, with a me- 
dian time interval of 14 months, and in the VSA data, time 
intervals range between 2 and 17 months, with a median 
time interval of 16 months. We have not quantified the vari- 
ability using Eqn. [2] as it is not suitable in the case where 
the level of variability is less than approximately the uncer- 
tainties on the flux density measurements and the number of 
data points is small; this is because it does not take account 
of the error on the mean flux density. 

Instead, we define the unbiased variability index for a 
pair of flux density measurements as follows: let Si and S% 
be the two measured flux densities with errors, including 
contributions from both calibration and thermal noise, ai 
and 02 respectively. We assume that, on average, the mea- 
sured change in flux density is related to the true change in 
flux density, ovar, by 



(Si — S2) = cr ins + a v 



(3) 



where crins = V^i+^J i s the measurement error on [Si — 
S2). The true change in flux density can therefore be ex- 
pressed as 



(4) 



We then define our unbiased variability index for a pair 
of flux density measurements as 



(•>) 



where S is the mean flux density. The factor 1/2S ensures 
that V lies between and 1, as is the case with V rms . V con- 
verges to Vims when the level of variability becomes large 



4 AMI Consortium: Franzen et al. 



compared with the uncertainties on the flux density mea- 
surements. 

When the value inside the square root in Eqn.[5]becomes 
negati ve, following a method similar to iBarvainis et al.l 
(2005), we set the unbiased variability index to be negative 
and this is given by 

V = -±y/\(S 1 -S2)*-W + <ri)\. (6) 
2.1.3 Results 

Figure [T] plots our unbiased variability index at 33 GHz 
against that at 16 GHz. It illustrates that the levels of vari- 
ability, for sources in this sample, are broadly correlated at 
the two frequencies. However, the plot cannot be used for 
detailed comparison of the variability at the two frequencies, 
because, for individual sources, the observations were carried 
out with different time intervals using the two telescopes. 

The median values of V at 16 and 33 GHz are 6.2 and 
4.7 per cent respectively. We stress that the median value 
of V at 33 GHz is unreliable because variability detected at 
this level with the VSA is not significant. Table [5] shows the 
distribution of V at 16 and 33 GHz for the sources in our 
sample. 

The most variable source from the AMI data is 
J1419+3822 with V = 64 per cent on a timescale of 524 
days, the source having decreased in flux density by about 
a factor of five; the corresponding value of V from the VSA 
data is 49 per cent. The most variable source from the VSA 
data is J 1727+4530 with V = 65 per cent on a timescale 
of 486 days, the source having increased in flux density by 
about a factor of five; the corresponding value of V from the 
AMI data is 51 per cent. 

We find a si milar level of variabi lity in our source sample 
to that found bv lSadler et al.l (|2006l ). This suggests that the 
two samples have very similar source properties, despite the 
difference of an order of magnitude in flux density levels. 

Given the level of variability that we find, over a signif- 
icantly shorter time period than separates our observations 
from those carried out by WMAP, it seems likely that the 
Eddington bias, found in the comparison of AMI/ VSA and 
WMAP flux densities discussed in Paper I, can resonably be 
attributed to variability. 

2.2 Source classification 

Following iBolton et ail (120061 ). we have used a x 2 -test to 
classify each of the sources as variable (VAR) or non-variable 
(NV) at both our observing frequencies. Our null hypothesis 
is that the source has a constant flux density and that the 
ith measurement of the source's flux density, Xi, is drawn 
from a Gaussian distribution with a mean, /i (which we take 
as the mean flux density measured for that source), and a 
standard deviation, ai, given by the uncertainty on the ith 
measurement of the flux density. The value of \ 2 is given, 
therefore, by 

X 2 = E^-^, (7) 

where n is number of observations made of the source. When 
the degree of confidence that the data do not support the 



null hypothesis is greater than 99 per cent we define the 
source as variable, otherwise we accept the null hypothesis 
and the source is defined as non-variable. We have included 
the classification for each source at both observing frequen- 
cies in Table [3] The mean separation between all pairs of 
observations for each of the sources is also included. We 
stress that, since this value varies from source to source, the 
timescale on which we sample variability also varies with 
source. 

At 16 GHz, of the 93 sources, for which we have more 
than a single observation, we classify 49 (53 per cent) as 
variable. At 33 GHz, we have more than one observation 
for 96 of the 97 sources in the sample; of these we clas- 
sify 30 (31 per cent) as variable. We note that if we use a 
90 per cent confidence level in the chi-squared test, the per- 
centages of sources classified as variable are 65 and 49 at 16 
and 33 GHz respectively. Of the 30 sources defined as vari- 
able at 33 GHz, 26 are also classified as variable at 16 GHz. 
This leaves four sources which are only classified as variable 
at 33 GHz; however, we believe that this is because the ob- 
servations of these four sources spanned somewhat different 
date ranges on the two telescopes (see Fig. 12. 2[) . The higher 
proportion of sources classified as variable at the lower fre- 
quency is, in large part, due to the larger errors on the flux 
densities measured by the VSA. We attempt, below, to make 
an assessment of the correlation in the variability between 
the two frequencies for those sources defined as variable at 
both 16 and 33 GHz, taking account of the measurement er- 
rors. In Fig. 12.21 we have provided the light curves of these 
sources. 

We have attempted to assess the variability as a func- 
tion of time, at both observing frequencies, for those sources 
which we have classified as variable. We have calculated the 
variability index (equations [5] and [6} for every pair of ob- 
servations for these sources. Figures [2] and [3] show plots of 
the variability index versus number of days separating the 
observation pairs at 16 and 33 GHz, respectively. The data 
have been averaged within bins of 50 days in width; the 
mean variability index is shown for observation pairs falling 
within each of the bins, along with the standard error on 
the mean. We have also carried out an identical procedure 
for those sources classified as non-variable; the results have 
been included in the figures for comparison. 

The plots indicate that, at both frequencies, the mean 
variability index for those sources defined as variable, in- 
creases significantly with time. In contrast, for the non- 
variable sources, there is little change in the mean variability 
index with time, which, in both cases, are close to zero. The 
mean variability index is generally higher for the variable 
sources at 33 GHz compared with that at 16 GHz. However, 
due to the larger measurement errors on the VSA observa- 
tions, only the most variable sources are selected as variable 
at 33 GHz; for the same reason, the level of variability re- 
quired for a source to be classified as variable is likely to 
be greater at 33 than 16 GHz. This could explain why the 
mean variability is generally greater at the higher observing 
frequency. 

Figure |4] shows the mean variability index versus mean 
flux density at 16 GHz, for sources that were observed mul- 
tiple times at that frequency. There appears to be no strong 
correlation between flux density and variability. We com- 
pared the mean variability index of the brightest 50 per cent 
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Figure 2. The mean variability index at 16 GHz versus number of 
days separating the observation pairs for those sources classified 
as variable from the AMI data (dashed) and non-variable from 
the AMI data (solid). The data have been averaged using bins of 
50 days in width. 
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Figure 3. The mean variability index at 33 GHz versus number of 
days separating the observation pairs for those sources classified 
as variable from the VSA data (dashed) and non-variable from 
the VSA data (solid). The data have been averaged using bins of 
50 days in width. 



of sources to the faintest 50 per cent; we found no statisti- 
cally significant difference in the mean variability index of 
the brighter sample compared to the fainter. Figure [5] shows 
a similar plot at 33 GHz. There are a number of fainter 
sources with low variability indices; this is, however, to be 
expected due to the significant contribution of thermal noise 
to the measurement uncertainties for the fainter sources at 
33 GHz. Again, we found no statistically significant differ- 
ence in the mean variability between the brighter and fainter 
sources. 




-5 5 10 15 20 25 30 35 40 45 
Mean variability index 

Figure 4. Mean variability index versus mean flux density at 
16 GHz. 



Mean variability index 

Figure 5. Mean variability index versus mean flux density at 
33 GHz. 



2.3 Correlation in the variability at 16 and 
33 GHz 

2.3.1 Methods 

We consider the group of 26 sources which are classified 
as variable from both the AMI and VSA data. For further 
analysis, we have required that a source has had at least two 
pairs of observations made simultaneously with the AMI and 
VSA; this produces a set of 16 sources. It was not possible to 
treat sources individually in this analysis, because, with only 
two or three pairs of simultaneous observations per source, 
the problem was too unconstrained. Instead, we treated the 
16 sources collectively, so that, with a total of 35 pairs of si- 
multaneous observations, the problem was well constrained. 
Figure [7] shows our joint dataset. Note that, before concate- 
nating the data, we normalized the fluxes so that the mean 
AMI and VSA fluxes of each source are equal to 1.0. The 
range of time intervals between observations is 51-277 days 
and the average time interval is 162 days. 

For this joint dataset, we then estimated the Pearson 
product-moment correlation coefficient R using a maximum- 
likelihood method. To carry this out, we form the likeli- 
hood of the data given covariance matrices for the signal 
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Figure 6. Flux densities at 16 (solid line) and 33 GHz (dashed line) of the sources which are classified as variable from the VSA data. 
Four of these sources (J0433+0521, J1058+0133, J1608+1029 and J2148+0657) are not classified as variable from the AMI data but this 
is thought to be because the observations of these four sources spanned somewhat different date ranges on the two telescopes. 
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Figure 6 - continued 
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and noise. Let Xi and j/i be the i th measurements of the flux 
densities at 16 and 33 GHz respectively, a x ,i and a Vt i the un- 
certainties on Xi and y% respectively and n the total number 
of measurements. The log-likelihood function is given by 

lnL(q) = (8) 
1 " 

constant — - [ln|Ci| + (x; — m)*Cj — m)] ,(9) 



where (x; — m)' is the transpose of (x^ — m), 
C, = S + N, 



Sn 
S21 



Sl2 
S22 



( Xi \ 


,m=( 


m x \ 


V Vi J 







the vector q 



' X,l 





Su 



(10) 
(11) 



(12) 

S12, S22) represents 



the parameters of a Gaussian model to be determined by 
maximizing the likelihood, the parameters m, and m y are 
mean flux densitiess at 16 and 33 GHz respectively, Su and 
S22 are variances of flux densities at 16 and 33 GHz respec- 
tively and Si 2 (= S21) is the covariance between flux densi- 
ties at 16 and 33 GHz. Then, R is a derived parameter and 
is given by 



R = 



Si 



V S11S22 

We used the nested sampling algorithm as i mple- 
mented in the MultiNest co de (|Feroz fc Hobsonl |200S| ; 
iFeroz. Hobson fc Bri dges 2008]) to sample points from ln(q) 
with the constraints — 1 ^ R ^ 1 and |d| ^ 0. 



2.3.2 Results 

We obtained a value of R — 0.955 ± 0.034 with a maximum 
likelihood value of 0.967. We repeated the calculation set- 
ting the noise covariance matrix term to zero. This yielded 
a value of R — 0.914 ± 0.030, i.e. instrumental noise leads 
to a significant decrease in the degree of correlation unless 
account is taken of the noise. This demonstrates the impor- 
tant role of the noise covariance matrix term in removing the 
bias in R resulting from instrumental noise. Figure [8] shows 
one-dimensional marginalized probability distributions for R 
with and without use of the noise covariance matrix term. 

Our analysis has shown that, for the variable sources 
in general, over time scales of several months, flux densities 
at 16 and 33 GHz are very highly correlated. Note that we 
have assumed that flux densities at 16 and 33 GHz are drawn 
from Gaussian distributions. This is unlikely to be the case 
for individual sources observed over long periods of time. 
However, in this particular case, we think it is a reasonable 
assumption to make given the fact that the analysis is run 
on a group of 16 sources with each source having approxi- 
mately the same number of observations. We have assumed 
that measurement errors at the two frequencies are Gaus- 
sian distributed and completely uncorrelated. The latter is 
certainly true, since we are dealing with two telescopes at 
different sites. Owing to the lack of data, it was not possi- 
ble to treat sources individually, making it hard to ensure 
that sources are given equal weight. We therefore stress that 
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Figure 7. Joint data-set, formed from pairs of simultaneous ob- 
servations at 16 GHz (solid line) and 33 GHz (dashed line) for 16 
variable sources, for which the Pearson product-moment correla- 
tion coefficient is estimated. 
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Figure 8. The one-dimensional marginalized probability distri- 
bution for the Pearson product-moment correlation coefficient 
with (solid line) and without (dashed line) use of the noise co- 
variance matrix term. 



our results are not necessarily applicable to any individual 
source but rather to the population. 

From this joint dataset, it is also possible to estimate 
the standard deviations of flux densities at 16 and 33 GHz, 
Si6 and S33 respectively, using this method. Again, these 
are derived parameters and are given by Sie = vftl and 
S33 = \fS£. We obtained a value of s w = 0.202 ± 0.028 
with a maximum likelihood value of 0.188 and a value of 
S33 = 0.224 ± 0.039 with a maximum likelihood value of 
0.198. Figure [9] shows one-dimensional marginalized proba- 
bility distributions for sie and s 33- 

These results suggest that there is little difference 
in the levels of variability at the two frequencies. Previ- 
ous studies suggest that there is only a slight increase in 
the degree of variability of compact radio source s with 
frequency. The ATCA was used by iTingav et all (|2003T l 
to observe a sample of 202 sources from the VSOP all- 
sky survey l|Hir abavash i et ail |2000| ) at 5 GHz over a pe- 
riod of 3-4 years. They measured a median variability 
at 1.4 GHz of 6 per cen t and at 8.6 GHz of 9 per cent. 
lOwen. Spaneer fc Cottonl l|l980h studied a group of sources 
with S90 GHz > 1 Jy and found only slightly more variability 
at 90 GHz than at 5 GHz. 
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Figure 9. The one-dimensional marginalized probability distri- 
bution for the standard deviations of flux densities at 16 (solid 
line) and 33 GHz (dashed line) 
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Figure 10. Mean spectral index versus mean variability index at 
16 GHz. 

2.4 Relationship between variability and spectral 
properties 

We measured spectral indices between 13.9 and 33.75 GHz 
using simultaneous observations with the AMI and VSA 
(Paper I). Figure [101 shows a plot of mean spectral index 
versus mean variability index at 16 GHz. As noted above, 
the timescale sampled by the variability index varies be- 
tween sources. Nevertheless, there is no reason to expect 
the timescale sampled for a source to be related to its 
spectral properties and, consequently, the plot is still use- 
ful for investigating the relationship between spectral in- 
dex and variability. The plot shows a trend for sources 
with more st e eply r ising spectra to be increasingly variable. 
iBolton et alj ([2006) have previously found evidence for a 
similar effect. They found no evidence for variability in a 
group of sources with a > 0.5, but 50 per cent of sources 
with a < 0.5 were found to be variable. 

We have attempted to quantify the correlation between 
spectral index and variability, in our data, by calculating the 
mean spectral index for the sources defined as variable and 
non- variable at 16 GHz. We find the mean spectral index 
for the variable sources as —0.06 ± 0.05 and for the non- 
variable as 0.13 ±0.04, where the uncertainties are standard 



errors on the means. Given the level of the errors, we find 
a significant difference between the mean spectral indices of 
the two populations. 



3 CONCLUSIONS 

In order to investigate the cm-wavelength- variability of flux 
density of sources selected at cm wavelengths, at the high- 
flux-density end of the source population, we have made ob- 
servations with the AMI (13.9-18.2 GHz) and VSA (33 GHz) 
of a complete sample of sources found with WMAP at 
33 GHz and with S33 ghz ^ 1.1 Jy. We applied a new al- 
gorithm to data obtained by near-simultaneous - typically 
a few minutes apart - observations at our two observing 
bands. This maximum-likelihood method allowed us to as- 
sess the correlation in the variability at the two frequencies, 
taking the measurement errors into account. We found that: 

(1) on timescales of « 1.5 years, 15 per cent of sources 
varied by more than 20 per cent at 16 GHz and 20 per cent 
varied by more than 20 per cent at 33 GHz; 

(2) this level of variability means that variability has to 
be taken into account in coping with source contamination 
in CMB work, particularly at £ > 2000; 

(3) this level of variability would indeed seem to explain 
the Eddington bias found in Paper I; 

(4) in a subset of sources classified as variable from 
both the AMI and VSA data, the Pearson product-moment 
correlation coefficient between the variability at 16 and 
33 GHz is 0.955 ± 0.034 and the degree of variability at 
16 GHz (0.202 ± 0.028) is very similar to that at 33 GHz 
(0.224 ±0.039); 

(5) there is a significant difference between the mean 
spectral indices of the variable sources (—0.06 ± 0.05) and 
non- variable sources (0.13 ± 0.04); 

(6) the AMI measurements show no strong evidence of a 
correlation between variability and flux density for the sam- 
ple as a whole. There appears also to be little difference be- 
tween the general l evels of variability in our sample and the 
ISadler et all (|2006l ) sample complete to 100 mjy at 20 GHz. 
These results seem somewhat surprising given that we do 
find a relation between spectral index and flux density. 



ACKNOWLEDGMENTS 

We are grateful to the staff of the Mullard Radio Astronomy 
Observatory for the maintenance and operation of the AMI 
and to the staff of the Teide Observatory, Mullard Radio As- 
tronomy Observatory and Jodrell Bank Observatory for as- 
sistance in the day-to-day operation of the VSA. We are also 
grateful to the University of Cambridge and PPARC/STFC 
for funding and supporting the AMI. We thank PPARC for 
funding and supporting the VSA project and the Spanish 
Ministry of Science and Technology for partial financial sup- 
port (project AYA2001-1657). MLD, TMOF, NHW, MO, 
CRG and TWS are grateful for support from PPARC/STFC 
studentships. 



10 AMI Consortium: Franzen et al. 



REFERENCES 

Akritas M. G., Bershady M. A., 1996, ApJ, 470, 706 
Barvainis R., Lehar J., Birkinshaw M., Falcke H., Blun- 

dell K. M., 2005, ApJ, 618, 108 
Bennett C. L., Bay M., Halpern M. et al., 2003, ApJ, 583, 

1 

Blandford R. D., Konigl A., 1979, ApJ, 232, 34 
Bock D. C.-J., Large M. I., Sadler E. M., 1999, AJ, 117, 
1578 

Bolton R. C, Chandler C. J., Cotter C, Pearson T. J., Poo- 
ley G. C, Readhead A. C. S., Riley J. M., Waldram E. M., 

2006, MNRAS, 370, 1556 

Cawthorne T. V., Rickett B. J., 1985, Nat, 315, 40 
Condon J. J., Cotton W. D., Greisen E. W., Yin Q. F., 

Perley R. A., Taylor G. B., Broderick J. J., 1998, AJ, 115, 

1693 

Davies M. L., Franzen T. M. O., Davis R. J. et al., 2009, 

MNRAS, in press 
de Zotti G., Ricci R., Mesa D., Silva L., Mazzotta P., Tof- 

folatti L., Gonzalez-Nuevo J., 2005, A&A, 431, 893 
Dennett-Thorpe J., de Bruyn A. G., 2002, Nat, 415, 57 
Dent W. A., 1965, Sci, 148, 1458 
Feroz F., Hobson M. P., 2008, MNRAS, 384, 449 
Feroz F., Hobson M. P., Bridges M., 2008, preprint 

(|arXiv:0809.3437|l 
Gonzalez-Nuevo J., Argiieso F., Lopez-Caniego M., Toffo- 

latti L., Sanz J. L., Vielva P., Herranz D., 2006, MNRAS, 

369, 1603 

Gregorini L., Ficarra A., Padrielli L., 1986, A&A, 168, 25 
Gregory P. C, Scott W. K., Douglas K., Condon J. J., 

1996, ApJS, 103, 427 
Griffith M. R., Wright A. E., 1993, AJ, 105, 1666 
Hinshaw G., Nolta M. R., Bennett C. L. at al., 2007, ApJS, 

170, 288 

Hirabayashi H., Fomalont E. B., Horiuchi S. et al., PASJ, 
2000, 52, 997 

Jauncey D. L., Kedziora-Chudczer L. L., Lovell J. E. J., 
Nicolson G. D., Perley R. A., Reynolds J. E., 
Tzioumis A. K., Wieringa M. H., 2000, in Hirabayashi H., 
Edwards P G., Murphy D. W., eds., Proceedings of the 
VSOP Symposium, Astrophysical Phenomena Revealed 
by Space VLBI, p. 147 
Lang K. R., Rickett B. J., 1970, Nat, 225, 528 
Linsky J. L., Rickett B. J., Redfield S., 2008, ApJ, 675, 413 
Lopez-Caniego M., Gonzalez-Nuevo J., Herranz D., Mas- 
sardi M., Sanz J. L., de Zotti G., Toffolatti L., Argiieso F., 

2007, ApJS, 170, 108 

Owen F. N., Spanger S. R., Cotton W. D., 1980, AJ, 85, 
351 

Rickett B. J., Coles W. A., Bourgois G., 1984, A&A, 134, 
390 

Rickett B. J., Lazio T. J. W., Ghigo F. D., 2006, ApJS, 
165, 439 

Sadler E. M., Ricci R., Ekers R. D. et al., 2006, MNRAS, 
371, 898 

Shapirovskaya N. Y., 1978, SvA, 22, 544 
Spangler S., Fanti R., Gregorini L., Padrielli L., 1989, A&A, 
209, 315 

Tingay S. J., Jauncey D. L., King E. A., Tzioumis A. K., 

Lovell J. E. J., Edwards P. G., 2003, PASJ, 55, 351 
Toffolatti L., Argiieso Gomez F., de Zotti G., Mazzei P. F., 



Franceschini A., Danese L., Burigana C, 1998, MNRAS, 
297, 117 

Watson R. A., Carreira P., Cleary K. et al., 2003, MNRAS, 
341, 1057 

Wright E. L., Chen X., Odegard N. et al., 2009, ApJS, 180, 
283 

Zensus J. A., Pearson T. J., 1987, Journal of the British 

Astronomical Association, 98, 48 
Zwart J. T. L., Barker R. W., Biddulph P. et al., 2008, 

MNRAS, 391, 1545 



Variability of WMAP point sources 11 



Table 3. Results for individual sources. The columns are the source name, taken from the NEWPS catalogue; the Equatorial 
coordinates (J2000), from the PMN or GB6 catalogues; the number of observations, n; the mean flux density, S; the source 
class, variable (VAR) or non- variable (NV); the mean variability index, V, with the mean number of days separating the 
observations shown in brackets and the spectral index, afl'g 5 . 
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0.80 


± 





35 


J1229+0203 


12 


29 


05.6 


+02 


03 


09 


3 


5 


29.20 


26.98 


NV 


NV 


-2.2 


(185) 


0.1 


(233) 


0.07 


± 





05 


J1230+1223 


12 


30 


48.8 


+12 


23 


36 


3 


5 


24.49 


17.67 


NV 


NV 


-2.8 


(184) 


0.1 


(233) 


0.41 


± 





05 


J1310+3220 


13 


10 


29.5 


+32 


20 


51 


4 


5 


1.21 


1.79 


VAR 


VAR 


13.4 


(228) 


12.8 


(212) 


-0.61 


± 





11 


J1331+3030 


13 


31 


08.0 


+30 


30 


35 


4 


5 


3.34 


2.06 


NV 


NV 


-2.7 


(228) 


-0.1 


(212) 


0.63 


± 





10 


J1347+1217 


13 


47 


33.4 


+12 


17 


17 


3 


4 


1.36 


1.03 


NV 


NV 


-2.8 


(270) 


4.1 


(252) 


0.62 


± 





20 


J1357+1919 


13 


57 


04.1 


+ 19 


19 


19 


3 


4 


1.98 


2.36 


NV 


NV 


1.5 


(292) 


1.4 


(251) 


-0.35 


± 





13 


J1419+3822 


14 


19 


45.9 


+38 


22 


01 


4 


4 


0.60 


0.86 


VAR 


VAR 


41.7 


(315) 


20.6 


(252) 


-0.08 


± 





18 


J1504+1029 


15 


04 


24.0 


+10 


29 


43 


3 


4 


1.60 


1.50 


VAR 


NV 


4.1 


(270) 


-1.2 


(251) 


0.00 


± 





17 


J1516+0014 


15 


16 


40.7 


+00 


14 


57 


3 


4 


1.01 


1.07 


NV 


NV 


-2.7 


(270) 


-4.8 


(252) 


0.06 


± 





22 


J1549+0237 


15 


49 


30.0 


+02 


37 


01 


3 


4 


2.16 


2.19 


NV 


NV 


-2.6 


(270) 


-7.4 


(251) 


0.08 


± 





15 


J1550+0527 


15 


50 


35.2 


+05 


27 


06 


3 


4 


2.82 


2.84 


NV 


NV 


1.9 


(270) 


1.4 


(244) 


0.02 


± 





10 
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Table 3 — continued 



J1608+1029 


16 


08 


46 


4 


+10 


29 


05 


3 


4 


1 


34 


1 


44 


NV 


VAR 





6 


(270) 


6 


8 


(244) 





11 


± 


0.19 


J1613+3412 


16 


13 


40 


8 


+34 


12 


41 


4 


4 


3 


21 


2 


68 


VAR 


NV 


6 





(315) 


5 


3 


(244) 





23 


± 


0.08 


J1635+3808 


16 


35 


15 


6 


+38 


08 


13 


4 


4 


2 


48 


2 


58 


VAR 


NV 


5 


9 


(315) 


2 


4 


(244) 


-0 


12 


± 


0.12 


J1638+5720 


16 


38 


13 





+57 


20 


29 


4 


4 


1 


92 


2 


45 


VAR 


NV 


2 


6 


(315) 


-1 


2 


(244) 


-0 


34 


± 


0.12 


J1642+3948 


16 


42 


58 





+39 


48 


42 


4 


4 


5 


30 


5 


26 


VAR 


VAR 


9 


5 


(315) 


9 


4 


(244) 





18 


± 


0.07 


J1651+0459 


16 


51 


09 


2 


+04 


59 


33 


3 


4 


2 


63 


1 


33 


NV 


NV 


-2 


3 


(270) 


-5 


3 


(244) 


1 


03 


± 


0.18 


J1720-0058 


17 


20 


29 


7 


-00 


58 


37 





5 




- 


3 


87 


- 


NV 








-3 


9 


(33) 






- 




J1727+4530 


17 


27 


28 


4 


+45 


30 


49 


4 


6 





57 





94 


VAR 


VAR 


31 


8 


(308) 


29 


8 


(237) 


-0 


67 


± 


0.17 


J1734+3857 


17 


34 


20 


5 


+38 


57 


45 


4 


4 





89 





93 


VAR 


VAR 


7 





(315) 


16 


3 


(244) 


-0 


21 


± 


0.17 


J1740+5211 


17 


40 


36 


6 


+52 


11 


47 


4 


4 


1 


02 





95 


VAR 


VAR 


15 


3 


(314) 


13 


5 


(244) 





12 


± 


0.17 


J1743-0350 


17 


43 


59 


2 


-03 


50 


06 


4 


5 


3 


80 


3 


66 


VAR 


VAR 


11 


1 


(245) 


11 


9 


(206) 





28 


± 


0.07 


J1751+0938 


17 


51 


32 


7 


+09 


38 


58 


3 


4 


5 


08 


6 


36 


VAR 


VAR 


33 


4 


(270) 


7 


2 


(244) 


-0 


19 


± 


0.07 


J1753+2847 


17 


53 


42 


5 


+28 


47 


58 


3 


4 


1 


78 


1 


85 


NV 


NV 


1 


7 


(293) 


-0 


9 


(244) 


-0 


14 


± 


0.13 


J1801+4404 


18 


01 


32 


2 


+44 


04 


09 


4 


4 


1 


40 


1 


47 


VAR 


NV 


6 


1 


(314) 


7 


5 


(244) 





04 


± 


0.15 


J1824+5650 


18 


24 


06 


8 


+56 


50 


59 


4 


5 


1 


28 


1 


34 


VAR 


NV 


11 


3 


(309) 


5 


5 


(226) 


-0 


26 


± 


0.12 


J1829+4844 


18 


29 


32 


1 


+48 


44 


46 


3 


4 


2 


70 


2 


32 


NV 


NV 


-2 


2 


(293) 


-5 


8 


(244) 





16 


± 


0.13 


J1955+5131 


19 


55 


42 


3 


+51 


31 


54 


4 


3 


1 


18 


1 


50 


VAR 


NV 


16 


5 


(314) 


-8 


4 


(227) 


-0 


63 


± 


0.14 


J1959+4034 


19 


59 


21 


8 


+40 


34 


28 





3 




- 


37 


43 


- 


NV 








-1 


8 


(319) 






- 




J2123+0535 


21 


23 


43 


4 


+05 


35 


14 


3 


3 


1 


41 


1 


41 


NV 


NV 


4 


9 


(278) 


-6 


3 


(229) 





00 


± 


0.19 


J2134-0153 


21 


34 


10 


4 


-01 


53 


25 


3 


2 


1 


90 


1 


75 


NV 


NV 


-2 


4 


(277) 


-3 


1 


(334) 





27 


± 


0.19 


J2136+0041 


21 


36 


38 


6 


+00 


41 


54 


3 


6 


5 


76 


4 


22 


NV 


NV 


-2 


6 


(183) 


1 


3 


(183) 





42 


± 


0.05 


J2139+1423 


21 


39 


01 


5 


+14 


23 


37 


3 


4 


2 


13 


1 


75 


NV 


NV 


3 


2 


(221) 


-7 


3 


(181) 





33 


± 


0.08 


J2143+1743 


21 


43 


35 


6 


+17 


43 


54 


4 


6 





50 





57 


VAR 


NV 


12 


9 


(206) 


-8 





(201) 


-0 


23 


± 


0.14 


J2148+0657 


21 


48 


05 


5 


+06 


57 


36 


3 


5 


5 


52 


5 


07 


NV 


VAR 


2 





(212) 


4 


3 


(206) 





19 


± 


0.05 


JZZUZ-j-^iZlO 


99 
ZZ 


09 
UZ 


1/1 


Q 
O 


1 AO 
+4Z 




oy 


A 
'1 


O 


O 


A 7 


A 
1 


oo 


VA Tt 


V-tVIX 


"1 9 
±Z 


n 
U 


^ZZo ) 




Q 
O 




— u 


9fi 
ZO 


□Z 


u.uo 


J2203+1725 


22 


03 


26 


7 


+17 


25 


42 


4 


5 


1 


18 


1 


25 


VAR 


VAR 


10 


3 


(228) 


7 


8 


(206) 


-0 


07 


± 


0.11 


J2203+3145 


22 


03 


15 


8 


+31 


45 


38 


4 


5 


2 


62 


2 


60 


NV 


NV 


-0 


7 


(228) 





7 


(206) 


-0 


08 


± 


0.08 


J2212+2355 


22 


12 


05 


9 


+23 


55 


31 


4 


4 





95 





96 


VAR 


NV 


3 


5 


(228) 


3 


6 


(179) 


-0 


14 


± 


0.12 


J2218-0335 


22 


18 


51 


8 


-03 


35 


40 


3 


4 


1 


65 


1 


33 


NV 


NV 


1 


8 


(213) 


-3 


7 


(181) 





47 


± 


0.15 


J2225+2118 


22 


25 


37 


6 


+21 


18 


17 


4 


4 


1 


11 


1 


27 


NV 


NV 





1 


(228) 


2 


3 


(179) 


-0 


07 


± 


0.11 


J2232+1143 


22 


32 


36 


6 


+11 


43 


54 


4 


4 


3 


79 


4 


52 


VAR 


VAR 


5 


7 


(252) 


3 


3 


(180) 


-0 


29 


± 


0.05 


J2236+2828 


22 


36 


20 


8 


+28 


28 


56 


5 


5 





99 


1 


60 


VAR 


VAR 


12 


6 


(234) 


8 


3 


(168) 


-0 


51 


± 


0.10 


J2253+1608 


22 


53 


58 





+ 16 


08 


53 


4 


5 


7 


17 


10 


73 


VAR 


VAR 


14 


8 


(228) 


22 


8 


(206) 


-0 


56 


± 


0.05 


J2327+0940 


23 


27 


33 


1 


+09 


40 


02 


4 


6 


1 


81 


2 


18 


VAR 


VAR 


6 


8 


(201) 


12 


8 


(201) 


-0 


27 


± 


0.08 


J2335-0131 


23 


35 


20 


1 


-01 


31 


14 


3 


6 





56 





79 


VAR 


NV 


15 


3 


(184) 


-0 


5 


(183) 





21 


± 


0.17 


J2354+4553 


23 


54 


21 


9 


+45 


53 


00 


4 


5 





98 





77 


NV 


NV 


-1 


7 


(228) 


-7 


4 


(206) 





32 


± 


0.14 



